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Abstract 

The eigen frequencies of radial pulsations of neutron stars are cal- 
culated in a strong magnetic field. At low densities we use the mag- 
netic BPS equation of state(EOS) similar to that obtained by Lai and 
Shapiro while at high densities the EOS obtained from the relativistic 
nuclear mean field theory is taken and extended to include strong mag- 
netic field. It is found that magnetised neutron stars support higher 
maximum mass where as the effect of magnetic field on radial stability 
for observed neutron star masses is minimal. 

*E-mail : vkg@ducos.ernet.in 
^E-mail : santokh@ducos.emet. in 
^E-mail : jda@ducos.ernet.in 
§E-mail : agoyal@ducos.emet. in 



1 



1 introduction 



It is well known that intense magnetic fields(B~10 12 ~ 13 G) exist 
on the surface of many neutron stars. Objects with even higher 
magnetic fields have been surmised and detected recently. Recent 
observational studies and several independent arguments link 
the class of soft 7-ray repeaters and perhaps certain anomalous 
X-ray pulsars with neutron stars having ultra strong magnetic 
fields, the so called magnetars. Kuoveliotou et al(1998) found a 
soft 7-ray repeater SGR 1806-20 with a period of 7.47 seconds 
and a spin down rate of 2.6 x 10~ 3 syr^ 1 from which they esti- 
mated the pulsar age to be about 1500 years and field strength 
of ~ 8x 10 14 G. Since the magnetic field in the core, according to 
some models, could be 10 3 — 10 5 times higher than its value on the 
surface, it is possible that ultra strong magnetic fields of order 
10 18 — 10 19 G or even 10 20 G exist in the core of certain neutron 
stars. The virial theorem, it is usually argued, gives an upper 
bound of about 10 18 G for the field inside a neutron star(Lai and 
Shapiro 1991). However it has also been claimed that this can- 
not be taken as the last word on the subject, since at the super 
high density inside the star core, general relativistic corrections 
to the virial theorem could increase the upper limit on the max- 
imum allowed magnetic fields substantially (Hong 1998). Ac- 
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cording to Kuoveliotou,a statistical analysis of the population 
of soft 7-ray repeaters indicates instead of being just isolated 
examples, as many as 10% of neutron stars could be magne- 
tars. It is therefore of interest to study the equation of state 
of nuclear matter and various properties of neutron stars under 
such magnetic fields. In this paper we study the radial oscilla- 
tions of neutron stars in the presence of super strong magnetic 
fields. Studies of radial oscillations is of interest since Cameron 
suggested (1965) more than three decades ago that vibrations 
of neutron stars could excite motions that can have interesting 
astrophysical applications. X-ray and 7-ray burst phenomenon 
are clearly explosive in nature. These explosive events probably 
perturb the associated neutron star and the resulting dynamical 
behaviour may eventually be deduced from such observations. 
Observations of quasi-periodic pulses of pulsars have also been 
associated with oscillations of underlying neutron stars (Chan- 
drasekhar 1964 a, b, Cutler et al 1990). 

The EOS is central to the calculation of most neutron star 
properties as it determines the mass range, the red shift as 
well as mass-radius relationship for these stars. Since neutron 
stars span a very wide range of densities, no one EOS is ade- 
quate to describe the properties of neutron stars. In the low 
density regions from the neutron drip density (~4xlO n ) and 
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upto /9 n ~3.0xl0 14 gm/cc the density at which the nuclie just be- 
gin to dissolve and merge together, the nuclear matter EOS is 
adequately described by the BPS model (Bayam, Pethick and 
Sutherland, 1971) which is based on the semi-empirical nuclear 
mass formula. We adopt this BPS EOS and its magnetised ver- 
sion as given by Lai and Shapiro(1991)in this density range. In 
the high density range above the neutron drip density p n , the 
physical properties of matter are still uncertain. Many models 
for the description of nuclear matter at such high densities have 
been proposed over the years. One of the most studied models 
is the relativistic nuclear mean field theory, in which the strong 
interactions among various particles involved are mediated by a 
scalar field a, an isoscalar-vector field u and an isovector-vector 
field p. Along with scalar self-intercation terms it can repro- 
duce the values of experimentally known quantities revelant to 
nuclear matter, viz. the binding energy per nucleon, the nuclear 
density at saturation, the asymmetry energythe effective mass 
and the bulk modulus and provides a good description Of nuclear 
matter for densities upto a few times the saturation density p c . 
In our study, we have used this nuclear mean field theory and 
its modification in the presence of a magnetic field. 
In section 2 a brief discussion of the EOS is given at zero tem- 
perature. In section 3 we present the formalism for pulsations of 
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the neutron star models computed here as a result of integra- 
tion of the relativistic equations. Section 4 deals with results 
and discussions. 

2 The Equation of State (EOS) for Nuclear 
Matter 

We shall describe nuclear matter at high densities by the rela- 
tivistic nuclear mean field model, including the /^-contribution, 
as extended to include strong interactions. For densities less 
than the neutron drip density we adopt the BPS model in the 
presence of magnetic field as developed by Lai and Shapiro. 

2.1 The Nuclear Mean Field EOS at High Densities 

We consider the charge neutral nuclear matter consisting of neu- 
trons, protons and electrons in /^-equilibrium in the presence of 
a magnetic field and at zero temperature (T = 0). The expres- 
sion for the various quantities like, number density, the scalar 
number density and energy-pressure etc in the [B = 0] case are 
well known. In the presence of magnetic field these expressions 
for charged particles are modified in a straight forward way, viz 

d 3 p = eB £ (2 - <U / dp, (1) 

spin u=Q 
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in the integrals appearing in the expressions for the various 
quantities listed above. The total pressure P and the mass en- 
ergy density (p) of the system are given by 

P = £ Pi + k—f < ^0 > 2 -k—r 2 M 2 - \bm n (g„af 

i=p,n,e ^ TTIlo ^ *J 

i=p,n,e ^ m uj ^ "la 

+ ^n(^) 3 + ^C(^) 4 + l -{^-f < po > 2 (3) 



In the above 

p„ = 



±- 2 [^y }n (2p% - 3m* 2 ) + m' n nn(^±p^)] (4) 



'p = 7"2 E (2 - MKp}, - mjtfn(!2-p)] (5) 

f4 "£'(2 - 5 »»)IW. - mlln(^^)\ (6) 



Pe = 



Pn = ^[P*nP}n( 2 Pn 2 ~ m n * 2 ) - m n *Hn( ^ n f — )] (7) 
ft, = ^ Z ( 2 " WMfp* + ^/V ^f 7 '* )) (8) 

p /-? v max II _|_ ft t 

= Tt E (2 - S„ )(PePfe + m e 2 /n(^±^)) (9) 

47T v= q Tfl e 
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m PiV * 2 = m p * 2 + 2ueB 
m e ,v = Trii + 2veB 



Pfp l^p ™pv 

*2 *2 *2 



Pfe 2 = ~ m 2 - 2ueB 
ml-m p = mn ~m n = -(— ) n s 
n s = n s p + n n s 



Z7T "Vt/ 

n B =n p + n n 

*3 



n 



3tt 2 



n P = 7T"2 S ( 2 - ^,o)P/p* 

p ^ Vmax 

n e = 7^2 E ( 2 - ^,o)P/e 
Tip — 7l e 

< A) >= ^OP - ™n) 

< cj >= n B 

fii= fH- ( ) < CJ > -J3i( — ) < A) > 



for i=n,p 



(27) 



For a given value of fi e and B we can find self consistently the 
values of m*, //*, //* and so the pressure and mass energy can 
be computed. This helps us to compute the equation of state in 
the mean field approximation as a function of fi e or ng. 

2.2 The Magnetic BPS model 

The total pressure of the hadron matter below the neutron drip 
is given by 



where £l is the bcc coulomb lattice energy and P e the pressure 
of the electrons in the presence of the magnetic field is given by 



P = P e (n e ) +P L = P e (n e ) + -e L (Z, n e ) 



(28) 



Pe = 



IIILLJl, . 

£ (2 - 8 v0 )Mel-m l e -2veB 




)] (29) 



(30) 



n e = —J £ (2 - ^o)k 2 -m 2 e - 2veBY 



(31) 
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Consider matter condensing into a perfect crystal lattice with a 
single nuclear species(A,Z) at the lattice sites. The energy density 
is 

e = ^W N (A, Z) + £ ' e (n e ) + e L (Z, n e ) (32) 

where Wn is the mass-energy of the nucleus(including the rest 
mass of nucleons and Z electrons) e e is the free electron energy 
including the rest mass of electrons. Following Lai and Shapiro 
we use for Wn, the experimental values for laboratory nuclei 
as tabulated by Wapastra and Bos(1976,1977). The elements 
taken in this paper are given below in the table along with their 
mass energy Wn(A,Z). At a given pressure P,the equilibrium 
values of A and Z are determined by minimising the Gibbs free 
energy per nucleon, 

e + P W N (A,Z) Z 2 4e L 

9 = — = — — + A { ^- m ' c) + ZM c (33) 

The neutron drip point is determined by the condition 

9mm = m n c 2 (34) 
Knowing A and Z the energy can be determined from Eq.(32) 



9 



3 Radial Pulsations of a Non-Rotating Neu- 
tron Star 



The equations governing infinitesimal radial pulsations of a non- 
rotating star in general relativity were given by Chandrasekhar(1964). 
The structure of the star in hydrostatic equilibrium is described 
by the Tolman-Openheimer- Volkoff equations 

dp -G(p + pc 2 )(m + ^) 



dr ^ cV(l-ff) 
dm 



(35) 



, = 4irr z p (36) 
dr 



dv 2Gm(l + ^ 



dr c 2 r(l — 



2Gm - 

(P'T / 



(37) 



Given an equation of state p(p), equations (35)- (37) can be in- 
tegrated numerically for a given central density to obtain the 
radius R and gravitational mass M = M(R) of the star. The 
metric used is given by 

ds 2 = -e 2 "c 2 dt 2 + e 2 V + r 2 (d9 2 + sin 2 9df) (38) 

If Ar is the radial displacement 

Ar 

? = — (39) 
r 

C = rVf (40) 
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and the time dependence of the harmonic oscillations is written 
as exp{iat) , one gets the equation governing radial adiabatic os- 
cillations(Chandrasekhar 1964,Datta etal 1998,Anand etal 2000) 



where 



^ F " -TT^ (42) 

G = -^[(W + 3,) + ^(r P )-?(r P )] (43) 

9 — — ^ e pyp + pc) o(-r) ( 44 ) 

j9 + /oc z r dr c 1 p + /9C Z ar 



A = -l„[l - (45) 
In the above equations T is the adiabatic index given by 

r = v -±^- d l (46) 

C^J9 dp 

The boundary conditions to solve equation (41) are 

C(r = 0) = 
£p(r = R)=0 (47) 

The expression for 5p as given by Chandrasekhar(1964) is 

5p(r) = - (48) 

G?r r 2 r 2 dr 

All these equations are totally model independent and are in- 

fact the same whether we are considering neutron stars, quark 
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stars or any other dense stellar object. The nature of the object 
being considered and the particular model affects the structure 
of the star and the frequency of radial pulsations only through 
the EOS. Notice that in Chandrasekhar (1964) and Datta etal 
(1998) the pulsation equations were written in terms of £ in- 
stead of C. Equation (42) along with the boundary conditions 
represent a Sturm-Liouville eigenvalue problem for a 2 . From 
the theory of such equations we have the well known results: (i) 
Eigenvalues a 2 are all real, and (ii) They form an infinite dis- 
creet sequence 

al<a\< o\ .... 
An important consequence of (ii) is that if the fundamental ra- 
dial mode of a star is stable (do 2 > 0), then all the radial modes 
are stable. 

4 Results and Discussions 

To study the structure and radial oscillations of neutron stars 
in the presence of a strong magnetic field we have employed the 
BPS model with its generalization in a magnetic field given by 
Lai and Shapiro(1991) below the neutron drip and the RMF 
theory above it. We have used the values of various couplings 
fixed by Ellis et al(1991) which provide the known values of 
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various nuclear matter parameters : 




0.003748 



0.01525Mev 



O.OllMev 



O.OllMei; 



-l 



-l 



-l 



c 



0.01328 



(49) 



As explained in section 2.1 and 2.2, for the RMF theory the 
equations are solved in a self-consistent manner for the effective 
masses and chemical potentials ,and hence the EOS. Below the 
neutron drip, the EOS is obtained by the minimization of Gibb's 
free energy as functions of A and Z. For this purpose we have 
employed 14 nuclei listed in the table. The problem is solved 
seprately for B = and 5^0. For each B this gives the EOS in 
the form P(ng) and p(ub)- The structure of the neutron star is 
then obtained from the integration of the oppenheimer-Volkoff 
equations. This then gives the profile of m,p and v as a func- 
tion of r for each star. One more quantity that is required is V 
which is then calculated directly from the EOS for all densities 
by using a quadratic difference formula for the derivative d ^ c 2y 
Along with the M-R relationship, one also obtains the gravita- 
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tional red shift 

Z = [1 - - 1 (50) 



c 2 r 



which can in principle be observed experimentally. 
The procedure to obtain the eigen frequencies is simple. We 
guess a value for a and integrate the equation outword from the 
centre upto the surface. The guessed value of a is varied till the 
boundary condition 

5p(r) = at r = R 

is satisfied. 

In figure 1. we plot mass in solar mass unit vs radius in Kms. 
for magnetic fields 0, 1x10 s , 5xl0 4 , 1x10 s MeV 2 (lMeV 2 = 
1.69 x 10 14 G) represented by the curves A, B, C and D respec- 
tively. It is worthwile to note that the magnetised neutron stars 
support higher masses. For very high magnetic fields the stars 
become relatively more compact. In figure 2. we present mass 
vs central energy density and in figure 3. we have plotted grav- 
itational red shift vs mass. In figure 4. a plot of time period 
of fundamental mode vs gravitational red shift is given for the 
same magnetic fields as in figure 1. It is interesting to note that 
for the observed neutron star mass (1.4 M ), the magnetic field 
has pratically no influence on radial stability. Similar trend is 
seen for the first excited mode as shown in figure 5. 
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TABLE: BPS Equilibrium Nuclei Below Neutron Drip 



Element 


BPS MASS-ENERGY 
(in units of 10 4 Mev) 


Fell 


5.2103 


mil 


5.7686 


mil 


5.9549 


Ni% 


6.1413 


Jxr 3Q 


8.0025 


O 84 

oe 34 


7.8170 


Gei 


7.6316 


Znf Q 


7.4466 


Nil! 


7.2621 


Ru\f 


11.7337 


Mo\f 


11.5495 


Zrlf 


11.3655 


S4? 


11.1818 


S4? 


11.3655 


Kr\f 


10.9985 
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Figure 1: Plot of mass in solar mass unit vs radius in Kms. for magnetic 
fields 0, lxlO 4 , 5xl0 4 , 1x10 s MeV 2 represented by the curves A, B, C and 
D respectively. 
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Figure 2: Plot of mass in solar mass unit vs energy density for magnetic 
fields 0, lxlO 4 , 5xl0 4 , 1x10 s MeV 2 represented by the curves A, B, C and 
D respectively. 
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Figure 3: Plot of gravitational red shift (Z) vs mass in solar mass unit for 
magnetic fields 0, lxlO 4 , 5xl0 4 , 1x10 s MeV 2 represented by the curves A, 
B, C and D respectively. 
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Figure 4: Plot of time period r for fundamental mode vs gravitational red 
shift(Z) for magnetic fields 0, lxlO 4 , 5xl0 4 , 1x10 s MeV 2 represented by 
the curves A, B, C and D respectively. 
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Figure 5: Plot of time period r for n=l mode vs gravitational red shift (Z) 
for magnetic fields 0, lxlO 4 , 5xl0 4 , lxlO 5 MeV 2 represented by the curves 
A, B, C and D respectively. 
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